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We use the density matrix renormalization group method to investigate the role of longitudinal 
quantized phonons on the Peierls transition in the spin-Peierls model. For both the XY and Heisen- 
berg spin-Peierls model we show that the staggered phonon order parameter scales as y/X (and the 
dimerized bond order scales as A) as A — > (where A is the electron-phonon interaction). This 
result is true for both linear and cyclic chains. Thus, we conclude that the Peierls transition occurs 
at A = in these models. Moreover, for the XY spin-Peierls model we show that the quantum 
predictions for the bond order follow the classical prediction as a function of inverse chain size for 
small A. We therefore conclude that the zero A phase transition is of the mean-field type. 

PACS numbers: 75.10.Jm, 71.38.-k 
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Since the 1930s the role of electron-phonon interactions 
on the stability of the metallic state in one-dimensional 
metals has been of abiding interest to both theoreti- 
cal chemists and physicists. In the 1950s Frohlichi and 
Peierls^ showed that as a consequence of electron-phonon 
interactions a one-dimensional metal is unstable with re- 
spect to an insulating (or Peierls) ground state. Equiv- 
alent conclusions were later derived for conjugated poly- 
mers by Ooshika£, and Longuet-Higgins and Salemi. All 
of these results were derived in the adiabatic limit (or 
Born-Oppenheimer approximation) in which the phonon 
degrees of freedom are regarded as the slow, classical 
variables. In this limit an infinite metallic system un- 
dergoes a Peierls transition to the lower symmetry insu- 
lating ground state for any value of the electron-phonon 
interaction. 

The question of whether such a conclusion remains 
valid for quantized phonons was first addressed by Frad- 
kin and HirscbA Using both renormalization group argu- 
ments and Monte Carlo simulations they suggested that 
quantum fluctuations do not destroy the Peierls state for 
non-interacting electrons. However, they also suggested 
that quantum fluctuations can destroy the Peierls state 
for non- interacting spinless fermions, giving a non-zero 
value for the critical electron-phonon interaction. 

More recent numerical calculations on various mod- 
els have also indicated that the Peierls state may be de- 
stroyed by quantized phonons. Caron and Moukour: 6 
used finite size scaling analysis of the spin gap to show 
that at a critical value of the electron-phonon interaction 
there is a Kosterlitz-Thouless transition in the XY spin- 
Peierls model with gapped, non- dispersive (or Einstein) 
phonons. Similar conclusions were made for the Heisen- 
berg spin-Peierls model, again with Einstein phononsi&. 
Finally, the Holstein model with Einstein phonons for 
both spinless 9 and spinfull 10 fermions, and the extended 
Hubbard-Pcicrls model with Einstein bonds phonons 11 
have also be shown to exhibit a Peierls phase transition 
at a non-zero value of the electron-phonon interaction. 

In this article we focus our attention on the role of 
quantized gapless, dispersive (or Debye) phonons on the 
Peierls state. We investigate both the XY spin-Peierls 



and Heisenberg spin-Peierls models. Spin-Peierls mod- 
els are the strong coupling limit of the Hubbard-Peierls 
model at half-filling. For open, linear chains the XY 
spin-Peierls model is equivalent to the non-interacting 
spinless fermion-Peierls model, allowing us to make di- 
rect comparisons to the essentially exact adiabatic pre- 
dictions. For both models we find that although quan- 
tized phonons reduce the amplitude of the broken sym- 
metry order parameter, in contrast to the previous in- 
vestigations, there is no evidence that the Peierls state is 
destroyed at a finite value of the electron-phonon inter- 
action. 

The Heisenberg spin-Peierls model with dispersive, 
gapless phonons is defined by, 



H = H 1 +H 2 
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ue and Pi are the conjugate displacement and momen- 
tum operators for the ^th site, and Se is the Pauli spin 
operator for that site. 

Defining the phonon frequency loq = y/K/M, the 
lattice degrees of freedom are quantized for dispersive 
phonons by introducing the phonon creation and annihi- 
lation operators, 
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respectively, where u = ^u)(jH. 
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TABLE I: The ground state energy (in units of J) of the 
Heisenberg spin-Peierls model as a function of the number of 
sites, TV. loq — J/v2 and g = 0.1. 









Number of oscillator levels per site 


N 


1 
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3 4 5 6 


16 


-7.142 


-7.330 


-7.458 -7.502 -7.518 -7.523 


28 


-12.437 


-13.031 


-13.286 -13.378 -13.415 -13.430 


40 


-17.746 


-18.733 


-19.116 -19.256 -19.314 -19.339 



TABLE II: The ground state energy, E, (in units of J) of the 
Heisenberg spin-Peierls model as a function of the density ma- 
trix eigenvalue product cutoff, e, the number of system block 
states, m, and the superblock Hilbert space size, SBHSS, for 
a 40-site chain with 5 oscillator levels per site, tuo = J/V2 
and g = 0.1. 
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We solve Eq. Q by the density matrix renormalization 
group (DMRG) method 10 i 12 i 13 i 14 for chains of up to 150 
sites. We use 5 oscillator levels per site, with typically 200 
block states and 100, 000 superblock states. One finite 
lattice sweep is performed at the target chain size. The 
convergence indicators are shown in Tables HI and ITU For 
linear chains with open boundary conditions we maintain 
constant chains lengths by fixing the position of the end 
sites. 

Wc first investigate the XY spin-Peierls model for arbi- 
trary phonon frequency and electron-phonon interaction 
on linear chains with open boundary conditions. For this 
model, 

S e ■ S e+1 -> - (S^S W + S^S^ +1 ) . (6) 

Using the Jordon-Wigner transformation the XY model 
on an open chain is equivalent to the non-interacting 
spinless model, where 

E x -> tJ2 (1 + - ( c W+i + H -C.) (7) 

e 

and t O- J/2. Following Fradkin and Hirsch£, we in- 
troduce the dimensionless electron-phonon interaction, 
g = y 'a 2 t/K = \n/2 (where A is the usual definition 
of the electron-phonon interaction^) . 

The Peierls broken-symmetry ground state is charac- 
terized by a non-zero value of the staggered phonon order 
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FIG. 1: The phonon order parameter versus the electron- 
phonon interaction for the spinless fermion- Peierls model. 
wo = t. The inset shows a linear extrapolation of Nm versus 
g to the origin. (The crosses indicate the evaluated points.) 
The calculations were performed on chains with fixed, open 
boundary conditions. 

parameter, defined by 

m(N) = ±-'£{Bt +1 -Bt)(-iy, (8) 

t 

where Bi = (b\ + bi)/2 is the dimensionless displacement 
of the £th site and N is the number of sites. To determine 
whether this order parameter vanishes as a function of 
the electron-phonon interaction in the asymptotic limit 
it is necessary to perform a finite size scaling analysis. 
Fradkin and Hirsch- suggested that the scaling relation 

m(N) = —F(Nm(oo)) (9) 
Ng 

should apply. Thus, curves of Nm(N) versus g are ex- 
pected to coincide at the critical value of g at which 
m(oo) vanishes. 

Figs and J2J show Nm versus g for the spinless 
fermion-Peierls model (Eqns QJ, © an d 0) f° r val- 
ues of loq — t and loq = lOt, respectively. It is clear in 
both cases that for small values of g m is linearly pro- 
portional to g for all chain lengths. Furthermore, the 
curves extrapolate linearly to the origin indicating a zero 
value for the critical g. Notice that since the distortion of 
the ith. bond from its uniform value, Sap, is proportional 
y / Uog(B£ + i — Be), the average staggered bond dimeriza- 
tion is oc g 2 or oc A for g <C 1. 

Since the spinless fermion model has an exact solu- 
tion the order parameter of the spinless fermion-Peierls 
model is readily obtainable in the adiabatic limit, where 
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FIG. 2: As for Fig. Q with too = Wt. 



the lattice displacements are treated classically. In order 
to ensure constant chain lengths the classical model is 
defined as, 



K 



(10) 



where T is self-consistently chosen so that J2e( u z ~ 
U£ + i) = 0. The ground state dimerization of Eq. i|l(J|) 
for linear chains is found iteratively using the Hellmann- 
Feynman procedure^. We define the normalized, stag- 
gered bond dimerization, Sq, as 
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where a is the average bond length. Fig. © shows So 
(scaled by g) of a 40-site chain versus g for the two quan- 
tum cases considered previously, as well as for the adia- 
batic limit. Evidently as g — > So/ g tends linearly — > 0. 
Thus, Sq oc g 2 oc A. Moreover, the phonon order param- 
eter decreases as a function of u>q. As g increases the 
deviation between the quantum and classical predictions 
increase. 

In the adiabatic and asymptotic limits the bond dimer- 
ization satisfies 



Sq = 4 exp 



1 



1 

71 A 



(12) 



for A <C 1. Here, n = 1 or 2 for spinlcss or spinfull 
fermions, respectively. However, there are significant fi- 
nite size corrections to this prediction as A — > 0. Fig. (JIJ 
shows the bond dimerization versus 1/N. Also shown 



FIG. 3: The scaled bond dimerization for a 40-site chain in 
the adiabatic limit (diamonds and full line), loo — t (circles 
and dotted line), and wo = Wt (triangles and dashed line). 
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FIG. 4: The scaled bond dimerization versus inverse chain 
length in the adiabatic limit. Crosses are the quantum result 
with luq = t. g — 0.01. 



are the quantum results for chains of up to 100 sites. 
For this small electron-phonon interaction there is excel- 
lent agreement between the quantum and classical pre- 
dictions. 

We now turn to a discussion of the Heisenberg spin- 
Peierls modeli&. Now the electron-phonon interaction 
g = y cn 2 J/K. Fig. 101 shows the phonon order param- 
eter versus g for ujo = J/y/2. The results are consistent 
with those described earlier, namely m cx g as g — > 0, in- 
dicating a zero critical value of g. (Qualitatively similar 
results are also found for ujq = lOJ/y/2.) 

The results described so far have been for linear chains 
with open boundary conditions. For these chains the 
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FIG. 5: The phonon order parameter versus the electron- 
phonon interaction for the Heisenberg-Peierls model. u>o = 
J/y/2. The inset shows a linear extrapolation of Nm versus 
g to the origin. 

ends break the energetic degeneracy between the other- 
wise equivalent 'A'and 'B' phases, which are related by 
a translation of one repeat unit. In order to check that 



open boundaries do not introduce any artefacts, we also 
calculate the phonon order parameter for cyclic chains. 
In this case it is necessary to explicitly break the de- 
generacy between the 'A' and 'B' phases by adding the 
symmetry breaking term 

p^S,.S, +1 (-l) £ (13) 

t 

to Eq. (1). Fig. JBJ) shows Nm versus p for a 100-sitc 
chain for various values of g. The results are remark- 
ably consistent with those for open chains. Extrapolat- 
ing Nm(p) to p = shows that Nm cx g, with values 
very close to those shown in Fig. (J5J. 

In conclusion, we have used the DMRG method 
to investigate the role of gapless, dispersive quantized 
phonons on the Peierls transition in the spin-Peierls 
model. For both the XY and Heisenberg spin-Peierls 
model we show that the staggered phonon order scales as 
g (and the dimerized bond order scales as g 2 ) as g — ► 0. 
This result is true for both linear and cyclic chains. Thus, 
we conclude that the Peierls transition occurs at g = 0. 
Moreover, for the XY spin-Peierls model we showed that 
the quantum predictions for the bond order follow the 
classical prediction as a function of 1/N for small g. We 
therefore conclude that the zero g phase transition is of 
the mean- field type, Eq. I|12l) . with a renormalized effec- 
tive value of A. 
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FIG. 6: The phonon order parameter versus the symmetry 
breaking parameter p for a 100-site cyclic chain, ujq = J/%/2. 



